Operational Time of Arrival in Quantum Phase Space 

Piotr Kochanski 

Centrum Fizyki Teoretycznej PAN, Al. Lotnikow 32/46, 02-668 Warszawa, Poland 

Krzysztof WodkiewiczP 
Instytut Fizyki Teoretycznej, Uniwersytet Warszawski, Hoza 69, 00-681 Warszawa, Poland 
0\ , (version time.tex; February 1, 2008) 

G\ 

^ ! Abstract 

An operational time of arrival is introduced using a realistic position and 
momentum measurement scheme. The phase space measurement involves the 
dynamics of a quantum particle probed by a measuring device. For such a 
measurement an operational positive operator valued measure in phase space 
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is introduced and investigated. In such an operational formalism a quantum 



£S) ' mechanical time operator is constructed and analyzed. A phase space time 

^^ ■ and energy uncertainty relation is derived. 
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I. INTRODUCTION 

The problem of the proper definition of the quantum mechanical time observable, and 
the physical interpretation of the associated uncertainty relation between time and energy, 
is a subject of a long lasting debate. This is due to the fact, that both in classical and 
quantum nonrelativistic mechanics, the time parameter is not an independent dynamical 
variable, and its definition requires the use of a nonlinear function of the canonical variables 
describing the system. In classical physics the form of this nonlinear function does not 
bring any conceptual difficulties, for instance, for a free particle of a unit mass and in one 
dimension, with an initial momentum p = p(t ) and an initial position q = q(to), the time 
of arrival from go to a fixed position q is given by 

t = t a + q -^i. (i) 

Po 

In the general case of an arbitrary motion it is necessary to invert the equations of motion 
for a particle q = q(t,q ,p ) in order to find the time parameter t = t(q,q ,p ). This 
function can be in general multivalued and the choice of the form t is selected according to 
the physical meaning of the given solution. 

In quantum mechanics the problem of the definition of the time observable associated 
with such a nonlinear transformation becomes much more complicated since canonical ob- 
servables are noncommuting operators. This makes, for instance, the definition like (|]) 
ambiguous due to the ordering problem of q and 4 operators. 

It was Pauli, who first originated the question whether one can define a time operator 
as a canonical variable conjugated to the energy of the system. In his book on quantum 
mechanics [T|, he pointed out, that although formally we can write [H, i\ = ih, this formula is 
unsatisfactory since the operator t cannot be Hermitian due to the semi - bounded character 
of the Hamiltonian spectrum. As a result, the uncertainty relation in the form AEAt > h 
has a disputable interpretation. 

One can derive, however, the above uncertainty relation using a Fourier decomposition 
of a nonstationary state in the form if)(x,t)=L dEi()(x,E) exp(—iEt/h). In such an ap- 
proach, the uncertainty relation is a simple consequence of Fourier analysis, similar to the 
one encountered in classical Fourier optics for time-dependent fields and their spectra. A 
similar argument can be applied to derive the position and the momentum uncertainty 
relation of a wave packet in wave mechanics. However, it is well known that this last un- 
certainty relation has a much deeper meaning being the consequence of the fundamental 
quantum incomplementarity of two canonical variables - momentum and position. Such an 
approach to quantum nonstationary states was presented first by Fock and Krylov 0, and 
in the basically same spirit by Mandelstamm and Tamm ||. In these works the relation 
between the energy of the system and its "inner" time |4j], which usually is the lifetime of 
the nonstationary quantum state, has been derived and discussed 

Another approach was proposed by Landau and Peierls and by Fock and Krylov 
who discussed the problem of the time-energy uncertainty relation analyzing a typical 
scattering experiment of a particle on a test particle. The essential conclusion of these 
contributions was that the uncertainty relation is connected with the inability of measuring 
precisely the energy of a given system in an arbitrary short period of time. 



Neither of the above proposals lead to the definition of the time operator, what's more, 
it has been shown by Aharonov and Bohm || that the interpretation given by Landau and 
Peierls is erroneous i.e., it is possible to measure the energy in a time period which violates 
the postulated uncertainty relation. In the same article Aharonov and Bohm proposed a 
different approach to the problem of the time observable. They argued that in order to 
measure the time of arrival to a certain fixed point one needs to have a "quantum clock" 
i.e., another quantum system, which reads the experimental data collected during a clock 
measurement. This assumption led to the first explicit definition of the time observable in 
the form proposed by Aharonov and Bohm: 

t = Uq\ + W) , (2) 

2 \ p p J 

where q,p are the operators of the clock particle. This definition formally leads to the 
commutation relation [H c , T c ] = ih, where H c is the Hamiltonian of a clock particle. However 
such a time operator has an important disadvantage, in an a priori imposed and physically 
unjustified choice of a symmetric ordering of momentum and position operators. On the 
other hand, such an operator has a number of interesting properties |||7| that indicate 
possible physical applications. 

This fact can be understood in two ways.. First, it has been shown in Ref. [§,0, that 
with the Aharonov - Bohm time operator (Q) one can associate a positive operator valued 
measure (POVM), or equivalently, a realistic measurement scheme, whose outcomes are 
described by such an observable. 

Another interpretation of the Aharonov - Bohm formula fl2|) might be given in the frame- 
work of the approach proposed by Kijowski ||10| . It was shown there, that by a construction 



of a probability distribution for a particle to pass at a certain moment a two dimensional ref- 
erence plane, one can define a time operator. The explicit formula for such a time operator 
is like d): 

Tk =sgn(p)- lq- + -q) , (3) 

2 V P P J 

where the multiplicative sign of the momentum operator assures the Hermitian character 



of the time observable [ll|. Physically, the sign function means that one distinguishes 
between particles moving to the left and to the right with respect to chosen two dimensional 
reference plane. For an initial beam of particles prepared in such a way, that it contains 
only particles moving either to the left or to the right, the time operator (^j reduces to the 
original Aharonov and Bohm formula @. One should however note that this time operator 
is built from the operators describing the measured particle, and not from some clock particle 
operators like in the Aharonov - Bohm approach. 

The time operator (||) is conjugated to the variable sgn(p)(p 2 /2), and as a consequence 



leads to a corresponding uncertainty relation [10,1] 



It has been shown in [12,13| that one can use the above definition as a starting point to 



look for a time observable. The definitions of the time observable given by Kijowski [10,1] 
and by Delgado and Muga [r2"|,|T3"|] are then equivalent. 



It shall also be stressed, that contrary to the claims of some authors [p|,p!2|, the time 
operator (§) is defined in the framework of standard quantum mechanics, although in a 



rather unusual representation of the wave function, which is, however, as good as any other 
representation (see comment [fll]). 

It is also worth noticing that a symmetrized operator like (0) is naturally connected with 
a probability distribution of arrival times derived from an expectation value of the properly 



defined positive current operator [p^pfj. 

After a number of early, classical papers, which stated the problem and revealed the main 
difficulties in dealing with the question of the time observable in quantum mechanics, there 
was a large number of further works. Although our brief review of the published literature 
is far from being complete, we would like briefly to describe several approaches in order to 
place our own work in some historical and logical context. An excellent study, of both older 
and newest researches connected with the problem of time observable in quantum mechanics, 
can be found in ]14] . 

Roughly speaking, the majority of the research efforts in this field can be divided into 
the following three approaches. 

In the first approach the main effort was concentrated on the proper definition and un- 
derstanding of various time and energy uncertainty relations. In such relations, the meaning 
of both time and energy, varied depending on the author and the presented method. The 
most complete, known to us, review of this problem can be found in Ref. [|15|. Newest 
investigations on this issue have been presented in Ref. 



16] 



Another way of dealing with the time observable was to define a time operator as a 
symmetrized function of the position and momentum operators |||7]]. This approach was 
discussed above and is justified by the work of Kijowski. 

The third method, the closest to our proposal, is based on the investigation of a particle, 
whose time of arrival to a given point is detected by a certain reference "quantum clock". 
This quantum clock is represented by another quantum system, which has some properties 
enabling us to read the time of arrival of the measured particle. Obviously, one can design 
many different models of such quantum clocks. Depending how one chooses the system 
representing the quantum clock, and how one models the interaction between the clock 
and the investigated particle different nonequivalent schemes of time measurement can be 
introduced. 



,16-18 



Examples of various "quantum clocks" can be found in 
Another interesting discussion of the realistic, irreversible, detector model has been pre- 

The decoherent histories approach to the arrival time problem has 



20|. In this paper a probability distribution to find 



sented by Halliwell [|TJ 
been given by the same author in Ref. 
a particle in a certain position in a given fixed time interval has been introduced, 
using similar idea presented earlier by Yamada and Takagi |21] . 



This has to be contrasted with the standard quantum approach, when such a probability 
is introduced only for a given moment of time. A similar approach, in which the detector is 
not specified as another quantum system, but rather is treated phenomenologically without 
explicit treatment of the interaction of the measured system and the experimental device, 



is given in Ref. [22 



One can also ask, if it is possible to choose among the different quantum clocks those 
which are the most precise, in a sense that they saturate the time - energy uncertainty 



relation. This question has been addressed in Ref. |23 



A very extensive and general analysis of various aspects of the proper treatment of the 
time observable has been given in the papers by Allcock |E3 and by Mielnik 



It is the purpose of this paper to define an arrival time operator based on a realistic 
momentum and position measurement. Our work is mainly motivated by a recent successful 
experimental approach, in which operational measurements of the quantum phase of opti- 
cal fields have been performed p6| |. Such operational measurements have been related to 



simultaneous joint measurements of position and momentum operators in a space extended 
by the measuring device. 

In this paper we define an operational time operator, which is connected with a very 
simple model of a joint momentum and position measurement. Such a measurement may be 
implemented |27|-50|, if additional degrees of freedom responsible for the measuring devices 



are involved. For such a measurement we introduce the operational time observable and a 
corresponding POVM. We show that for a special case of the quantum clock state, the arrival 
time operator involves an antinormal ordering of the creation and annihilation operators 
forming the canonical variables. The operational phase space measurement leads to a time 
and energy phase space. For such a phase space and for an operational time of arrival a 
time and energy uncertainty relation is derived. 

This paper is organized in the following way. A short review of the general operational 
approach will be discussed in Section II. In Section III we present a quantum model for a 
joint measurement of position and momentum in a space extended by the measuring device. 
This allows to construct a corresponding POVM relevant to the definition of the operational 
time observable. From this POVM a positively defined operational probability distribution 
function called propensity is derived. This propensity leads to a mean time of arrival. In 
Section IV we introduce an explicit form of the operational time of arrival. We derive the 
time operator and show that the specific ordering of the position and momentum operators, 
forming this observable, depends on the properties of the measuring device. Using the time 
operator, constructed in such a way, we derive in Section V an operational time and energy 
uncertainty principle. In Section VI we introduce a time and energy phase space. We show 
that an arrival time distribution can be interpreted in terms of a positive and negative flux 
flow. Finally Section VII contains some concluding remarks. 

Through the rest of the paper we shall work in one dimension and use the units in which 
h — 1, m — 1. 



II. OPERATIONAL QUANTUM MEASUREMENT 
A. Operational POVM 

In order to describe realistic experiments involving measuring devices it is sometimes 



necessary to go beyond the idealized measurement scheme proposed by von Neumann |3T 
In this approach it is postulated, that statistical outcomes of a measurement of a certain 
observable A with eigenvalues and eigenvectors A\a) = a\a) are described by the spectral 
measure: 

^(a) = |(a|^)| 2 , (4) 

where \tp) G 7i is the state vector of the measured system. It is known that the spectral 
measure contains all the relevant statistical informations about the investigated system but 



it makes no reference to the apparatus employed in the actual measurement. Due to this 
property A will be called an intrinsic quantum observable. In order to provide a more 
realistic approach to the quantum measurement we need to analyze carefully the dynamics 
of the combined system involving both the measured system and the measuring apparatus. 
We shall call for short the measuring apparatus and all its associated devices a filter device. 



Such an approach has been presented by many authors |27H30| and we shall use here the 
formulation and notation from |p7| , p8|] . 

Let the state of the measured system be described by the density operator g in a Hilbert 
space Ti and the state of filter device by gjr in a Hilbert space Tijr. The evolution of the 
combined system, including the interaction between the filter and the investigated system 
necessary for the measurement to happen, is given by a unitary operator U(t, to) acting in 

n®n T : 

Q(t ) <g> y{t ) — ► U(t, t )g{t ) <g> gr{t )U\t, t ) , (t > t ). (5) 

In order to probe the system state with such an interaction, the statistical readouts of 
the filter device are collected. The corresponding probability distribution corresponds to a 
propensity of the measured state to take on certain states of the filter. Due to the interaction 
the propensity of the probed state to be in one of the filter states \a) is 

Pr(o) = Tr sys {gTrrig^\a)(a\U}} = Tr sys {pf(a)}, (6) 

where J 7 (a) is a filter dependent POVM satisfying the normalization condition 

dafi(a) = 1. (7) 



B. Operational operators 

In view of the linear relation between the propensity and the POVM, the operational 
statistical moments of the measured quantity 

daa n Pi(a) = (AP) (8) 

define uniquely a set of operational operators 

AP = I daa n F(a). (9) 



A real number a is a classical quantity which enables us to read out the result of the 
quantum mechanical measurement. Since a real measurement operates in a classical world 



the outcomes of such an experiment are known only if they are macroscopically recorded |32 
This is why the result of such an experiment will be given by the probability distribution of 
a certain classical variable. In the presented approach this distribution is just the quantum 
propensity Pr(a). 



The operators @ are called operational operators, because they represent quantities 
measured in a real experiment involving a dynamical coupling of the measured system with 
the filter. We see that in a realistic measurement, with a filter, the spectral decomposition 
of an observable A is effectively replaced by a POVM daJ r (a) [p9fl . As a rule, the algebraic 
properties of the set A^ differ significantly from those of the powers of A. For instance, 
(Ay) 2 ^ A\p ■ This property will have important consequences in the discussion of the 
uncertainty introduced by the measurement. As an example the operational spread of a 
measured quantity, described by the of operational operators, is given by 

5 A = (AP) - (ig } > 2 = H 5 - a 2 . (10) 

This relation will play an important role in the formulation of the operational uncertainty 
relation. 

It is seen that the propensity and the operational operators are natural extension of the 
spectral probability distribution and the intrinsic operators. The difference is that the oper- 
ational observables carry information about the system under investigation and the selected 
measuring device. There always exists a physical mechanism (a realizable experimental pro- 
cedure) generating any desired POVM. This is guaranteed by the Naimark extension, and 



the reduction of the projection measure into the Hilbert space of the measured system p9 
According to this theorem one can always extend the Hilbert space with the defined POVM 
to a space with a projective measure. 

III. JOINT MOMENTUM AND POSITION MEASUREMENT 

A. Joint measurements in phase space 

In this section we present a simple model of a joint position and momentum measurement 
proposed originally by von Neumann |31j and elaborated later in Ref. fl33|-|35|j . 



The interaction Hamiltonian between the system (a particle) and two filter particles 
(labeled by 1 and 2) is given by 

H I = 5(t)(qp 1 -pq 2 ). (11) 

After the pulsed interaction with the system, the measured readouts q_M and pm of the filter 
variables q\ and f>2 are: 

<zi — > qm = U^qiU = q + qi- (72/2 , 

V2 — >Pm = U^p 2 U = p + p2 -pi/2 . (12) 

In the combined space of the two filter particles, we introduce filter operators q^r = g 2 /2 — qi 
and fir = f>2 — Pi/2. In the product space of the system and the filter the two measured 
observables (|12"D are: 

QM = q-qT and ]3.m = j3 + p r . (13) 

We recognize in the relations (|T^) the position and the momentum observables of a composed 
system introduced and discussed first by Einstein, Podolsky and Rosen in their famous paper 
on the completeness of quantum mechanics 



In our case the filter particles provide an example of the Naimark extension of q and p. 
In the extended space, the filter readout variables commute: 

[qm,pm] = [q-qt,p + pA = °, ( 14 ) 

and, as a result, can be measured simultaneously. From this definition, we can introduce in 
this extended space a time observable: 

f M = f^ , (15) 

Pm 

where the ordering of the commuting measured observables is irrelevant. This expression 
will be central in the following investigations. The reduction of this operator to the Hilbert 
space Ti of the measured particle will provide an operational definition of the time of arrival 
in phase space. As an example we obtain in the space of the particle the following operator: 

T = TrJM. (16) 

Vpm) 

Using the definition of the propensity (|6]) from Section II (with a being q and p), and the 
above interaction Hamiltonian we obtain that the propensity is a phase space distribution: 



Pr(?,p) = 7T 



dxexp (—ipx)ijj(x)J-'*(q + x) 



(17) 



where ip(x) and ^{x) are the wave functions of the particle and the filter in the position 
representation labeled by x. In fact the filter wave function is an overlap of the wave 
functions of the two particles 1 and 2 forming the measuring device. The POVM in the 
Hilbert space of the system that corresponds to this measurement scheme is 

F(Qi p) — — ex P (i<lP — WO) \J : ) (F\ ex P (~iqp + ipq) ■ (18) 



B. Mean time of arrival 

We shall apply this measurement to probe the time of arrival of a freely moving particle 
probed by a filter device. In this case the mean time of arrival will be simply related to 
the mean relative position of the measured particle with respect to a fixed filter position. 
In order to describe such a measurement we shall apply the derived above phase space 
distribution. 

In such a case the propensity becomes time dependent via the wave function tp(x,t) of 
the measured system. In order to carry all the calculations in an analytic form, we shall 
assume that the time measurement is performed on a freely moving particle described by 
the following Gaussian wave function: 

/^^ /(OT*-*°)n (19) 



characterized by the initial mean position Xq, the initial mean momentum k$ and width S. 
The filter wave function, at the measurement time (t = 0), is assumed to be in a stationary 
state with a Gaussian profile centered around go with width a 



2 



■no 2 



1/1 



2/^2 



F{x) = -— e -^-t«r/°\ (20) 



As it is seen, due to the pulsed interaction between the filter and the particle, the dynamics 
of the filter particle is unimportant. 

For this choice of the two wave packets the time evolution of the propensity is easily 
calculated 



. 1 vW ( ^a^ + a^ip-kof 

Pi{q,P,t) = - ====== == exp 



n ^/At 2 + (5 2 + a 2 ) 2 V W + (S 2 + ° 2 ) 2 

( -2 (5 2 (q -q + x + k t) 2 + a 2 (q-q + x + ptf) \ 
X 6XP V At 2 + (5 2 + a 2 ) 2 ) ■ [ZL) 

To understand the meaning of the propensity associated with this measurement we calculate 
its first two phase space moments 



q(t) = q ~ (x + k t) , (22a) 



, t 2 5 2 + a 2 

~ l -i ■ 

S 2 



^) 2 =?+-> + —T- . (22b) 



P = k , (22c) 



1 1 
6 2+ V 2 



P 2 = k 2 + T2 +- 2 - (22d) 



We see, for instance, that the average q{t) gives the relative position of a particle mea- 
sured with respect to the filter position. This fact reflects the fundamental property of the 
measurement in which the arrival time to a fixed position is monitored by a filter device. 

The cross section of the propensity in q for a fixed value of p, as a function of t is 
depicted on Fig. [TJ. In this case the particle moves from a point x towards the detection 
point go of the filter. Note a full symmetry between positive and negative values of t, around 
the measurement time set at t = 0. With the passage of time, the distance between the 
measured particle and the filter particle is shrinking, and finally it is zero, meaning that the 
measurement corresponding to the arrival time at the point go has been performed. Due 
to the unavoidable dispersions of the particle and the filter wave packets the propensity 
is spreading in time, deteriorating the measurement precision. The accuracy of such a 
measurement might be improved in several ways. We can use well prepared wave packets 
of the system i.e., wave packets with a small dispersion, or we may prepare the filter state 
in a squeezed state [J3fi|], with a reduced position dispersion. The measured momentum is 
constant in time, which is a unique feature of the Hamiltonian with a pulsed (Dirac delta 
function) interaction in time. 



IV. TIME OBSERVABLE 

A. Operational time of arrival 

From the phase space propensity derived in the previous section it is straightforward to 
define a mean arrival time and its moments using the following definition: 



(p) = J r dqdp {l) Pr(9 ' P ' t} ' (23) 

where the integration region T is the whole phase space. According to the general rela- 
tion between the moments @ and the operational operators ©, we introduce a family of 
operational time operators: 



(f ( n) ) =(-) (24) 



where in terms of the operational time POVM (18): 



fM=Jdqdp(£\ F(q,p). (25) 

We notice immediately that for our choice of the propensity the p-integral is divergent due 
to a singularity at p = 0. This singularity can be easily understood on the physical ground, 
related to the concept of the arrival time discussed above. The physical reason for the 
singularity is that the particle with zero momentum will never arrive to the measurement 
point, implying an infinite arrival time. 

Obviously if we change the measurement scheme in such a way that (introducing for 
example a proper external potential) the particle is forced to move towards the detection 
point, the infinity in (|25|) would be removed. In fact such a divergence is not unusual in 
quantum mechanics, when for instance, vectors of the Hilbert space corresponding to an 
infinite mean value of the position operator are used. For the time operator operationally 
based on the measurement of the arrival time, we shall select only such states for which the 
probability of finding the particle with zero momentum is very small. 

In order to remove the divergence we can apply a regularization procedure, removing 
from the integration region the values of p near zero. As a result we obtain: 

7f) = J dqdpOj) F(q,p), (26) 

where for e > T £ = {(q,p) \ q G M.,p e R\ [— e, e]}. A similar regularization for a differently 
ordered time operator was described in great detail in Ref . || . Our conclusion is analogous 
to the results of this paper. If the initial wave packet is prepared in such a way that it has a 
vanishing distribution of momentum near p = 0, the final results involving the time operator 
are independent of the regularization procedure. The condition which must be fulfilled is 
that \k \5 > 1. 

On Fig. H we have depicted the phase space propensity for a selected set of values of 
t,q ,x ,5,a and k . From this plot of the propensity, we see that it is enough to set, for 

10 



example, ko = 40 and 5 = 0.1 (in our units) to have a propensity practically equal to zero 
in the required region. It is clear that if both the particle wave function and the filter wave 
function are badly localized, the notion of the time of arrival looses its meaning. In this case 
a different time operator ought to be introduced. 

Because the only significant contribution in the phase space integration comes from the 
momenta concentrated around p = ko we can expand 1/p in (^6]) in the series around ko, 
(from a mathematical point of view such an expansion of the integrand will lead usually, 
after integration, to an asymptotic series). This greatly simplifies the integration in (p6|), 
and as a result we obtain the following leading contributions to the two first moments: 

(f «> = (*Z-^ - t) + 0(l/(ky), l/(ky), l/{kl5a)) , (27a) 

(f ( 2 )> = (f «) 2 + i^L + o(i/(ky), i/(ky), i/(a») , (27b) 

As we see, the leading contribution reproduces the semiclassical approximation of the motion 
of the particle. The physical interpretation of the obtained results is clear. The mean of, 
the introduced time operator, gives the time of flight between the two points x + k t 
and go- This time of arrival is given as a ratio of the distance between two points and 
the particle momentum. Obviously there are quantum corrections to the formulas (p7|). 
These corrections are much more pronounced when a second moment of the time operator 
is investigated. 



B. Ordering and time of arrival 

We shall derive now a more direct formula for the operational time operator (|25| ) in the 
case of a particular choice of the filter wave packet parameters. We express the position and 
the momentum operators, in terms of the oscillator creation a) and annihilation a operators, 
using the standard relations q = (l/\/2)(a) + a) and p = (i/\/2)(aJ — a). We set the filter 
wave function to be a state corresponding to (q\a\J-') = 0. Such a state is a Gaussian wave 
function with dispersion a = \/2 in our units). Using this filter state and these operators 
we can rewrite the POVM as follows 

F{a) = -D{a)\0)(0\D\a) = \a)(a\ with fd 2 aF(a) = i, (28) 

where D(oc) = exp (aa) — a* a) is the Glauber displacement operator in the complex a = 
(l/v / 2)(g + ip) phase space |]57|] . In such case, i.e., for such a filter wave function, the 
propensity is simply P(a) = ~(a\p\a). We recognize in this expression the so called Q - 
representation used in quantum optics applications. It is well known that this function is 
related to the antinormal ordering of the creation and annihilation operators (a) to the right 
and a to the left) [|7j. 

The time observables then take a very simple form 

(a\ n 
l J ^ (29) 

11 



where the triple dots denote the antinormal order of the creation and annihilation operators. 
We see that the measurement scheme with this particular choice of the filter wavepacket 
leads naturally to an antinormal ordering of the creation and annihilation operators, which 
reflects a specific ordering of the momentum and position operators forming the operational 
arrival time observables (0). For a different choice of the width parameter a, we will have 



a different ordering of these operators (see [34] for details) . 



As it was discussed above, only these states, which have vanishing amplitudes for mo- 
menta equal to zero are physically interesting, and for such states the time operator written 
above is well defined. Such states specify the domain of the operator (PD|). 

In that spirit we will find a more explicit formula for the first two moments of the time 
operator. We need to express the inverse of antinormally ordered (in the discussed above 
sense) momentum operator as a direct function of momentum operator centered around 
k ^0: 

.1. 1 . 1 . 1 v^, 1N „. fp-k 



!-£(-!)"= FP) :■ (») 



In this formula k will be understood as mean value of momentum of a detected particle, 
and that states with zero momentum are excluded. Naturally, the above equality holds only 
when we assure that there is a domain of vectors T> = {ip G 7i | ||(p//eo — 1)^11 < 1}- All 
the following formulas will have a physical meaning for states belonging to the this domain. 
In the next derivations we shall use the following algebraic identity: 

: exp(Ap): = exp ( — + Xp J , (31) 

where A is an arbitrary parameter. Using this formula simple algebra gives 

[ f = ¥o T,{^) H„«P-M), (32) 

where H n (z) are Hermite polynomials. In order to find a second moment of a time operator 
we need the following expression 



p 2 e^o de p - e kl ^-J \ 2ik ) 



(n+l)H n (i(p-A;o))- (33) 



All higher powers :(l/p) n : are easily obtained by the subsequent derivation over e, and then 
taking the limit e — > 0, as shown above for n — 2. 

Combining these algebraic properties we obtain the following formulas for the first two 
moments of the operational time operator: 

\ lb \J fi \J J 



12 



where :1/j3 2 : was calculated above. 

In the subspace D an arbitrary moment of the operational time operator can be written 
in the following form: 

oo 

f^ = J2 f M (36) 



fc=0 



where the subscript index is of order ((p—ko)/ko)) +1 , and a systematic method of calculating 
all Tr fc ™ follows from the algebra presented above. For sufficiently large values of ko i.e., for a 
fast moving particle we can well approximate the formula for the time operator by the first 
few terms in the series. For T^ we have: 



,(i) , ^(i) , ^(i) , /n I fp-k 



f« = f,' + f' + r 2 ' + o r-i^ ,A , (37) 
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[1] T [2] T "\\ h J 'A* 
where 

T$ = I , (38a) 



f{i) _ 1 ( q P-k p-kp q 

[1] -"2U^T ^TF 1 ' (38b) 



T M = 2 (^ (—J + (—J S ' ' <38C) 

The zero order term reproduces the classical limit corresponding to an arrival time of a 
particle with momentum fc at a fixed position q. 
For T^ we have 

?2 



,(2) _ g _,_ i_ 



T M ~ M + 2H ' (39a) 



Tff = -A (^ + #^ - ^ (V^ + ^V) • 09b) 



[1] A; 2 V k ^ H k V 2A; 2 V/ ^o ^ fc 

Even in the zero order approximation the second moment of the time observable does not 
reproduce the classical limit, in particular we see that (T^) 2 ^ T^ 2 \ 

Usually when we deal only with mean values of observables, the difference between results 
obtained with the help of an intrinsic operator (if such exists) and the operational operator 
is rather simple and straightforward. The first moment of the operational operator equals 
to the intrinsic operator, or differs by a term simply connected to the properties of the of 
the filter device. The situation becomes quite different for the second moments. In the 
operational formalism the second moment contains a contribution representing quantum 
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fluctuations of the filter. Due to the additive character of the presented measurement the 
time operator T^ contains an additive quantum noise contribution from the filter system. 
We conclude this section investigating the eigenfunctions of the time operator. The 
eigenvalue problem for the operational time of arrival has the following form: 

f {l) \Xr) = r\x T ). (40) 

In the following we derive the eigenvectors for several approximated values of T^\ In the 
zero approximation, for TL , the eigenstate can be simply calculated in the momentum 
representation. For a given eigenvalue r this eigenvector is just 

x (°)(p) = (p\ T ) = e - iTko{p - ko) = e~ iTk o , (41) 

where we have introduced a new variable T = r{p — fc ) Ao- As it is seen, the eigenstates 
are labelled by the momentum of the inspected particle. If we treat r as a physical time, 
we can write rp = q, rk = q and then T can be interpreted as the time of flight from the 
point q to a fixed point q. The spectrum of this operator is continuous. 

The eigenfunctions of the time operator involving the higher order contribution TL can 
be obtained solving a first order differential equation. As a result we obtain: 



exp ( iTk log f 1 



p-fco 



X?> (P) = ' V " r ^V ±^~ • (42) 

This eigenfunction will belong to the domain T>, if (p — ko)/ko < 1. The connection with 
the eigenfunction in the lower order of approximation and the interpretation of this result 
is best seen upon expanding the eigenfunction in powers of the parameter (p — ko)/ko, we 
have then 

^ )K ( 1 + ^ + ...) exp (-^( 1 + i^ + i(^) 2 + ...)),( 4 3, 

where, as above, T = r(p — k )/k . From this expression we conclude that the phase of this 
eigenfunction is the classical arrival time T with nonclassical contributions resulting from 
the higher order term TL . 

V. OPERATIONAL TIME-ENERGY UNCERTAINTY 

In the approach presented in this paper it is easy to give a clear meaning to the time and 
energy uncertainty relation. This is because that, on the operational level, we can describe a 
measurement of both the time of arrival for a given system and its energy. The operational 



time operator is given by Eq. (25), and the corresponding operational moments of energy 
of the system are: 

EW = Jdp (^j F(p), (44) 
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where F(p) = J dqF(q,p) is the momentum marginal of the phase space POVM (|18|) . It 
should be pointed out that the defined below uncertainty relation is given between the 
measured time of arrival and the measured energy of the chosen system, and as a result, the 



criticism of Aharonov and Bohm Q] and Peres [30| concerning those approaches where one 
tried to construct the uncertainty principle between the time of arrival detected by the filter 
and the energy of the investigated particle (two commuting quantities) does not apply here. 
Using the general formula for the operational spread introduced in Section II, we intro- 
duce the operational spread of the time operator 

5T = (fW) - (T«> 2 . (45) 

This quantity can be calculated using Eq.( |27b| ). The corresponding operational spread of 
the energy is: 

5E = (E®) - (E^) 2 . (46) 

The mean operational energy (& 1 ') has been already calculated fl22d| ) and we need only 

'5 2 + a 2 \ f5 2 + a 2 \ 2 



<*">-s[«+*(t£)*»(t£)' 

Simple algebra gives, as a result, an operational time-energy uncertainty in the form: 



(47) 



5E 5T = ((E^) - ((E^) 2 )((fW) - (f «) 2 ) > -^-tf_ . (48) 

2 da 

When both 5 and a are approaching simultaneously zero, the right hand side of the above 
inequality tends to 1 (i.e., h), meaning that the "conventional" uncertainty relation is re- 
produced if the particle and filter states are very well prepared. However, we have to 
remember that these calculations have been performed under the assumptions that k 5 ^> 1 
and k a ^> 1. This corresponds to the position spread of the wavepackets very small and 
hence giving the momentum spread very large. 

We conclude this section calculating the commutator involving the energy & 1 ' and the 
operational time operator T^ derived for a filter state leading to an antinormal form. As 
it was said before, these first moments of the operational operators should in some sense be 
the closest to the intrinsic time and energy observables. A simple algebra shows that: 

|f< 1 »,^»] = ! ( 1 + H_^);_l E¥ ;. (49) 

The right hand side of this commutator is complicated, but for states belonging to the 
domain T> we obtain that the leading term is: 

[f^,E^]=i + o(t^l\. (50) 

Only in this limit, one can associate the traditional interpretation of T^ and E^ as canon- 
ically conjugated variables for time and energy. 
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VI. TIME AND ENERGY PHASE SPACE 

A. Probability distribution of time of arrival 

The definition of the operational time of arrival operator is related to the following 
moments (q/p) n of the corresponding phase space propensity. It seems natural then to 
introduce a differently parameterized propensity i.e., instead of working with Pv(q,p,t) as 
a function of q, p we shall define a new arrival "time" variable 6 = q/p (to be not confused 
with the time t labelling the evolution of the system wave function): 



Pr (9, t) =s(--e\= J dqdp s(- - 6] Pr(q, p, t) . (51) 

From the properties of the Dirac delta function we obtain 

Pr(0,t)=/ dppPr(6p,p,t)+ dppPi(-6p,-p,t) =Pr+(0,t) H- Pr _(0, t) . (52) 
Jo Jo 

The interpretation of this result is simple in terms of the right (+) and left (— ) operational 
probability flux. 

This result should be contrasted with the standard quantum current operator 

3' 0*0 = 2 (P\ x )( x \ + \ X )( X \P) ■ (53) 

The expectation value of this current seems to be a natural and intuitive candidate for the 
probability distribution of the time of arrival at the given space point x. It is however 
difficult to make a practical use of this intuition since the flux ([53]) is nonpositive and it is 
difficult to associate with such a nonpositive flux a probability distribution. A discussion of 



this problem has been presented by Delgado [13 



Because the position and the momentum of the particle cannot be specified with ar- 
bitrary precision there is a possibility for a quantum backflow contribution to (j(x)), if it 
is understood as a probability distribution. As it was pointed out in Ref. |L3||, when this 
backflow part of the probability distribution becomes negligible one can use this formula 
to define the time observable. Indeed, if a wavepacket is constructed in such way that the 
mean value of the momentum is positive and large, and a small negative part of the flux 
can be neglected, the probability distribution for an arrival time is approximately positive. 
However, this backflow contribution is interesting in itself E3] since it is of purely quantum 



mechanical origin. This is why it is tempting not to exclude it from the flux considerations, 
and in our approach based on the operational propensity, there is no need to do so. 

Contrary to the quantum mechanical flux (|53|), the operational flux given by ([52]) is 
clearly positive. And as we have seen in the previous sections an association between this 
flux and the time observable is very natural and simple. In fact the flux ( |5"2]) is exactly the 
corresponding time of arrival probability distribution, 

(T (1) ) = f d66Pi{6,t). (54) 
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Obviously, in our approach we deal with the two parts of the flux. The "positive" part 
Pr + (#, t) corresponds to the particle moving in the direction given by the mean value of the 
momentum k . The "negative" (or backflow) part has obviously the opposite meaning. The 
presence of these two parts in the flux is clearly a purely quantum feature, because in the 
classical world the particle moves either to one or to the other direction with probability 
one. For our choice of the wavepackets the negative and the positive parts of the propensity 
(E2I) might be written explicitly: 





Pr±(M 


1 V<? 2 $' 

2vr 

x exp < - 


V4t 2 + (a 2 + 5 2 ) 2 

A(0,t) 

-2k 2 (6 2 T cl (t) + a 2 ' 

At 2 + (a 2 


(2+^V(erf(e) 

r d (0) + A(0,0))\ 
+ 5 2 ) 2 j' 


-! 


where 


















T d (t) = 


_ qo- 


(x + k t) 
k 


T d (0) = q °- Xo 
k 


) 








A(6,t) = 


--o\e 


+ t) 2 + 6 2 f 


+ U 2 a 2 {5 2 + a 2 ) 

A 


, A(0,0) = 


-5 2 a 2 (5 2 
4 v 


+ a 2 




f = 


2k ( 


j\Q + t)T d (0) + 5 2 6T cl (t) + A(0, 0)) 







^/2A(6,t)^4t 2 + {a 2 + 5 2 ) 



(55) 

(56) 
(57) 
(58) 



and erf(z) is the error function. When the mean momentum k$ is large and positive the 
contribution of the negative flux Pr_(#, i) becomes small. In this case the ratio of these two 
parts might me approximated as follows: 

Pr+(M) 2V5FU \e)j l ' 

When the particle is moving fast, k ^> 0, then £ becomes large and the above ratio tends 
to zero. Finally, the formula for the complete flux is given by 



1 v^ 


2 5 2 ^At 2 + (a< 


+ 5 2 ) 2 


2tt 


A(9,t) 




1 \G 


2 6 2 ^/4t 2 + (a r < 


+ 5 2 ) 2 



Pr( ^ t} " ^^ 6XP < At 2 + (a 2 + 52)2 > ~ 



, crf(0c::p f -^ 2 (^-^(t)) 2 ] 
2vr A(6,t) ^ er HUexp| 2A( ^ t) |. [W) 

This quantity is depicted on Fig. |^ for a particular choice of the wave function parameters. 
As it's seen the backflow part of the flux is much smaller then the dominating "positive" 
part. What is important is the fact that in our considerations we do not have to neglect 
the "negative" momenta contribution. In fact our model measurement scheme allows to 
measure this part of the probability distribution. 

It is also interesting to note, that having defined the time observable Delgado and Muga 



l4|13ll , also have managed to associate with it a positively defined current. This shows that 
such a relation is quite universal. 

In analogy to the time distribution, we can define a complementary quantity, a probability 
distribution connected with the energy measurement: 

17 



Pr(E) = s(^-E\=J dqdpS (^- - E) Pr(q,p,t) . 



(61) 



Due to the pulsed interaction of the filter with the system, this distribution is time inde- 
pendent. Indeed, a simple calculation shows that: 

Pr(E) l ,pglexpl ^ 2 (kl + E) } ( sV(k V2E) \ 



B. Time and energy phase space 

In conclusion of this section we shall introduce a combined time and energy phase space 
for the operational measurement. In such a phase space we have a joint energy and momen- 
tum distribution defined as: 



Pr(0, £) = *(£- 1?W| -*) . (63) 

The time distribution ( (52|) and the energy distribution (|6T|) are marginals of this joint time 
and energy distribution. For the Gaussian wave functions we calculate 



pi{e E) = 2 vW cogh I 4k V2E(5 2 T d (t)6 + a 2 T d (0)(9 + t) + A(0, 0)) 



7T y/At 2 + (a 2 + 5 2 ) 2 \ At 2 + (a 2 + 5 2 ) 2 

2E 



"2fc 
x exp ' 



At 2 + (a 2 + 5 2 ) 2 



T A(0, t) + 5 2 T cl (t) + a 2 T d (0) + A(0, 0) 



'0 



(64) 



The meaning of this quantity becomes clear, when we look on its plot depicted on Fig. [|. We 
see the characteristic separation into two parts, one corresponding to the positive momenta 
(relative to ko) and a much smaller contribution from the negative part. Further, we notice 
that for small values of energy the propensity is narrow in the E direction but rather broad 
in the 9 direction. This shows that the time measurement looses its meaning for very slow 
particles (small ko). The situation changes in a complementary way when the energy is 
increased. This is depicted on Fig. ^ where the propensity becomes very sharp in the 9 
direction. This assures a sharp measurement of the time of arrival. In this case the energy 
measurement becomes less accurate. 

The fact, that we have these two limits of the propensity is of course a manifestation of 
the time-energy uncertainty principle. 

VII. CONCLUSIONS 

With the help of the operational formulation we have introduced an operational time 
observable associated with a specific measurement scheme. Such an operational approach 
is not universal, in a sense, that the specific form of the time observable depends on the 
quantum state of the filter device. This operational approach allows for a natural and clear 



definition of the time of arrival. Our result is an intuitive quantum counterpart of the classi- 
cal time of flight measurement between two points. Based on such an operational approach 
the time - energy uncertainty principle has been introduced. The idealized measurement 
scheme allows to provide a link between the quantum mechanical flux of the operational 
propensity and the time observable. This has been done without neglecting the backflow 
part, which is very interesting. 

The time operator discussed in this paper is related to a specific measuring scheme. 
The question remains what have we learned about the intrinsic time operator from our 
discussion. Our view is that various fundamental time operators will have many features of 
the operator T^ . In fact we have shown that in terms of the commutation relations, classical 
limit and the physical interpretation, this operator has many features of an intrinsic time 
observable. However the problem is that we cannot provide a first principle derivation of this 
observable from some fundamental assumptions. The operational observable has been build 
using a specific detection procedure, but its overall properties should describe a intrinsic 
time operator in a reasonable way. 
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FIG. 1. Plot of the evolution of the cross section of the propensity Pi(q,p = 40) for momentum 
p = 40 at different times, with go = 3.5, xq = 1, 5 = 0.1, a = 0.1, ko = 40 in units h = m = 1. 




FIG. 2. Plot of the propensity Pr(g,p) at the time t = 0.01, with go = 3.5, xq = 1, S = 0.1, 
a = 0.1, ko = 40 in units h = m = 1. Near p = 0, the propensity is practically zero justifying our 
approximation . 
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FIG. 3. Plot of the propensity Pr(0,t) as a function of 9 at t = —0.1, with qo = 3.5, xo = 1, 
5 = 0.1, a = 0.1, ko = 10 in units h = m = 1. The right hand part of the plot shows the "positive" 
momenta contribution, whereas on the left hand side part of the plot we see the contribution from 
the "negative" momenta — the backflow effect. 




FIG. 4. Plot of joint propensity Pr(E, 9) at the time t = 0.01, with q$ = 3.5, xq = 1, <5 = 0.1, 
a = 0.1, &o = 40 in units H = m = 1. Again we see that there are two parts, since there are 
contributions from both "negative" and "positive" momenta. 
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FIG. 5. Plot of joint propensity Pr(E, 9) at the time t = 0.01, with go = 3.5, xq = 1, 5 = 0.1, 
a = 0.1, ko = 40 in units h = m = 1. The range of that plot is chosen such that we see only 
the positive momenta contribution. For low energies the time of arrival is very inaccurate — the 
propensity is very broad, but the measurement of energy is precise; for large momenta the situation 
is reversed, which is the consequence of the uncertainty relation for time and energy. 
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